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Improved  Correlations  for  Augmentor  Static  Stability 

Darin  A.  Knaus*,  Patrick  J.  Magari1,  Roger  W.  Hill',  and  Scott  D.  Phillips® 

Creare  Incorporated,  Hanover,  NH  03755 

and 

Barry  V.  Kiel  " 

Air  Force  Research  Laboratory,  Wright-Patterson  Air  Force  Base,  OH  45433-7607 

Flame  stability  is  critical  to  the  operational  performance  of  combustion  systems  in 
propulsion  and  power  generation.  Current  predictive  tools  for  flame  stability  are  based  on 
deeades-old  empirical  correlations  that  have  limited  applicability  for  modern  combustor 
designs.  Recent  advances  in  computational  fluid  dynamics  (CFD)  and  advanced  combustion 
diagnostics  have  provided  new  insight  into  the  fundamental  processes  that  occur  in  these 
flows.  Reaeting-flow  CFD  has  yet  to  mature  to  a  level  where  it  can  be  practically  applied  as 
a  design  tool  for  this  problem.  This  paper  describes  a  new  methodology  for  analyzing  flame 
stability  intended  to  provide  designers  with  a  significantly  improved  near-term  predictive 
capability.  Our  predictive  methodology  is  based  on  a  Damkohler  number  ( Da )  approach. 
Simplified  CFD  calculations  are  used  to  calculate  relevant  flow  timescales,  and  reactor 
model  calculations  are  used  to  characterize  the  important  chemical  timescales  in  the  system. 
These  timescales  are  used  to  form  a  Da  number  that  is  used  to  determine  stability.  Unlike 
existing  stability  correlations,  which  are  typically  based  on  global  parameters,  the  new 
methodology  predicts  stability  based  on  the  local  Da  field.  In  addition,  the  new  methodology 
is  expected  to  reduce  the  requirement  for  on  empirical  inputs  when  developing  a  stability 
model.  The  methodology  is  demonstrated  using  a  simple  rig  configuration.  Future  work 
will  involve  application  of  the  methodology  to  more  complex,  realistic  augmentor  conditions. 


Nomenclature 


CFD 

=  Computational  Fluid  Dynamics 

D 

=  Flame  Holder  Characteristic  Dimension 

Da 

=  Damkohler  Number 

£ 

=  Turbulent  Dissipation 

k 

=  Turbulent  Kinetic  Energy 

LBO 

=  Lean  Blow  Out 

LES 

=  Large-Eddy  Simulation 

P 

=  Combustor  Inlet  Pressure 

PDF 

=  Probability  Density  Function 

PSR 

=  Perfectly-Stirred  Reactor 

PaSR 

=  Partially-Stirred  Reactor 

RANS 

=  Reynolds-Averaged  Navier  Stokes 

Re 

=  Reynold’s  Number 

St 

=  Strouhal  Number 

Tad 

=  Adiabatic  Flame  Temperature 

^ chem 

=  Chemical  Timescale 

Tflow 

=  Flow  Timescale 

T 

1  o 

=  Initial  Temperature  for  Ignition  Delay 
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U  =  Combustor  Inlet  Velocity 

u  ’  =  Turbulence  Intensity 

®  =  Equivalence  Ratio 

I.  Introduction 

Flame  stability  is  critical  to  the  operational  performance  of  combustion  systems  in  propulsion  and  power 
generation.  Existing  design  systems  in  industry  today  are  based  on  parametric  testing  and  generic, 
simplified  correlations.  Empirical  correlations  for  predicting  flame  stability  parameters  such  as  lean  blowout  (LBO) 
are  often  based  on  global  inputs,  such  as  inlet  temperature  and  mean  flow.  This  type  of  approach  understandably 
has  limited  accuracy,  especially  for  state-of-the-art  systems  that  push  the  performance  envelope. 

Recent  advances  in  Computational  Fluid  Dynamics  (CFD)  and  combustion  diagnostics  have  provided  new 
insight  into  the  significant  processes  affecting  augmentor  stability.  Recent  data  have  allowed  flame  dynamics  near 
blowout  to  be  directly  observed  with  high  spatial  and  temporal  resolution.  For  example,  the  interaction  of  vortices 
shed  from  an  augmentor  flameholder  with  the  recirculation  zone  has  been  identified  as  a  key  process  associated  with 
blowout  and  both  static  and  dynamic  stability  under  certain  conditions 1’2,3,4.  CFD  has  provided  new  insight  into 
these  flows;  however,  the  technique  is  often  impractical  to  apply  as  a  design  tool  due  to  computational  cost  and 
complexity. 

The  objective  of  this  research  program  is  to  develop  computationally  tractable  reduced-order  models  for 
predicting  augmentor  static  stability  that  incorporate  insight  into  physical  processes  obtained  via  advanced 
diagnostics  and  CFD  calculations.  These  models  will  provide  combustor  and  augmentor  designers  with  a  new  tool 
that  can  help  bridge  the  gap  between  outdated  correlations  based  on  global  inputs  and  complex  reacting  flow  CFD 
calculations  and/or  expensive  parametric  experimentation. 

II.  Background 

The  design  of  combustion  systems  existing  in  industry  today  is  based  on  empirical  testing  and  generic,  simplified 
correlations.  Current  tools  have  limited  applicability  to  modern  combustor  designs6.  Several  of  the  correlations  in 
use  are  analyzed  in  detail  below.  Experimental  data  for  LBO  fuel-to-air  ratio  are  collected  and  a  correlation  for 
predicting  LBO  is  developed  based  on  the  global  parameters  varied  in  the  experiments,  such  as  inlet  velocity.  Many 
critical  factors  are  not  incorporated  into  the  correlations,  and  the  correlations  generally  do  not  extrapolate  well 
beyond  the  parameter  space  on  which  they  were  defined.  Not  surprisingly,  these  correlations  often  poorly  predict 
LBO,  especially  when  used  to  predict  state-of-the-art  augmentor  designs  that  push  the  performance  envelope  and 
utilize  unique  geometries. 

The  cost  of  full-scale  qualification  of  combustion  systems  has  increased  substantially  such  that  parametric 
testing  of  design  iterations  to  solve  operability  and  performance  shortfalls,  has  created  significant  system 
development  costs.  A  design  system  capable  of  predicting  the  influence  of  the  geometry  and  design  features  on 
operability  and  performance  could  save  millions  of  dollars  in  program  development  costs  and  help  ensure  that 
engine  systems  are  highly  optimized.  An  improved  design  system  must  capture  the  influence  of  dimensionality  and 
local  interactions  that  drive  the  behavior  of  the  system. 

CFD  techniques,  such  as  RANS  and  Large-Eddy  Simulation  (LES),  have  been  applied  to  solve  reacting  flows  for 
simplified  configurations7"13.  A  recent  AIAA  special  issue  was  dedicated  to  this  topic5.  CFD  methods  for  reacting 
flows  have  significant  computational  cost,  so  that  they  are  not  currently  a  practical  tool  for  designers.  However, 
reacting  flow  CFD  can  be  used  to  generate  a  large  quantity  of  spatially  and  temporally  resolved  data  that  cannot  be 
readily  obtained  experimentally.  For  this  reason,  reacting  flow  CFD  is  useful  for  developing  insight  into  these  flows 
that  can  be  applied  to  reduced-order  models  that  are  computationally  tractable. 

The  purpose  of  this  research  program  is  to  bridge  the  gap  that  exists  between  current  design  tools  and  available 
diagnostic  and  computational  capabilities.  Our  approach  involves  use  of  tractable,  mature  computational  techniques 
to  calculate  several  timescales  relevant  to  stability.  The  timescales  are  used  to  form  a  Damkohler  number,  Da,  that 
is  used  to  predict  stability.  Damkohler  number  approaches,  based  on  global  inputs,  have  been  developed  in  the  past. 
The  primary  distinction  between  previous  efforts  and  the  current  research  is  that  local  Da  values  in  the  regions  that 
determine  stability  are  calculated.  Below,  we  review  existing  stability  correlations  present  an  improved 
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methodology  for  predicting  stability,  demonstrate  the  methodology  using  LBO  data  for  an  augmentor  rig,  and 
discuss  challenges  and  next  steps  for  this  approach. 

III.  Existing  Stability  Correlations 

When  early  researchers  began  studying  these  flows  in  the  1950’s,  they  primarily  approached  the  problem 
through  direct  experimentation14'16.  Researchers  conducted  parametric  studies  based  on  global  parameters,  such  as 
inlet  velocity  and  temperature.  They  then  used  their  experimental  results  to  derive  stability  correlations  for  the 
prediction  of  blowout  equivalence  ratio  as  a  function  of  global  inlet  parameters.  A  loose  physical  justification  for 
the  correlation  was  typically  given;  however,  they  are  highly  empirical  in  nature.  Because  experiments  are  costly 
and  time  consuming,  the  parameter  space  on  which  early  stability  correlations  are  based  is  limited.  We  refer  to  these 
correlations  here  as  the  empirical  class  of  correlations. 

Two  commonly  used  empirical  correlations  are  those  of  DeZubay14  and  King16.  The  basis  for  DeZubay’s 
correlation  was  a  series  of  rig  experiments  using  disk  shaped  flameholders  in  a  circular  duct.  DeZubay’s  rig  was 
subscale  with  a  diameter  of  7  cm.  DeZubay’s  experiments  all  involved  hopane  fuel.  Inlet  velocity,  pressure,  and  the 
diameter  of  the  disk  were  varied.  Velocity  was  varied  from  12  to  168  m/s  and  pressure  was  varied  from  0.2  to  1  atm. 
Based  on  the  measure  LBO  data,  DeZubay  derived  a  correlation  parameter  for  predicting  LBO  ®.  DeZubay’s 
parameter  is  of  the  form: 


p  0.95  0.85  v  ’ 

More  recently  other  researchers  have  modified  the  DeZubay  parameter  to  include  a  term  for  inlet  temperature. 
When  plotted  against  this  parameter,  DeZubay’s  LBO  data  collapsed  to  a  curve,  a  stability  loop,  that  can  be  used  to 
predict  lean  and  rich  blowout. 

King’s  empirical  correlation  is  based  on  conditions  that  are  more  representative  of  an  actual  augmentor.  JP-4 
fuel  was  used  and  measurements  were  made  at  elevated  inlet  temperature.  King’s  experiments  spanned  inlet 
velocities  from  122  to  198  m/s,  pressures  from  0.35  to  0.85  atm,  and  inlet  temperatures  from  700  to  1,033  K.  King’s 
rig  was  full-scale  with  a  diameter  of  65  cm.  A  single  flameholder  geometry  was  considered,  so  that  no  geometry 
information  appears  in  the  correlation.  King’s  correlation  parameter  is: 

pO.324rl.O7(750_  ^  0.252  (2) 

King’s  correlation  predicts  LBO  only,  and  the  stability  curve  has  a  much  different  shape  from  DeZubay’s  curve  due 
to  the  way  that  the  parameter  is  formed:  DeZubay’s  parameter  increases  with  U,  whereas  King’s  parameter 
decreases  with  increasing  U. 

The  limitations  of  the  empirical  correlations  are  numerous.  The  correlations  are  derived  based  on  a  limited 
parameter  space  and  do  not  generally  apply  outside  of  that  parameter  space.  For  example,  DeZubay’s  correlation 
cannot  be  expected  to  perform  well  for  non-disk  shaped  flameholders  or  configurations  with  widely  varying 
blockage  ratios.  King’s  correlation  clearly  cannot  predict  a  configuration  with  U  >  750  (fps  units  are  used  when 
evaluating  the  correlation).  Neither  correlation  can  account  for  different  fuel  types  or  other  potentially  significant 
parameters  such  as  inlet  turbulence  intensity.  These  types  of  correlations  are  used  by  designers  as  a  starting  point 
but  cannot  be  expected  to  accurately  predict  new  configurations.  Modified  versions  of  the  empirical  correlation 
parameters  have  been  developed  by  engine  manufacturers  but  the  modified  forms  are  considered  proprietary  and  are 
not  openly  disclosed. 

More  recently,  researchers  have  improved  the  physical  descriptions  used  to  develop  stability  correlations  with 
the  hope  of  developing  more  accurate,  general  stability  correlations.  In  the  improved  physical  description,  stability 
is  seen  as  being  determined  by  competition  between  a  flow  and  chemical  timescale  (rc/,em).  We  refer  to  this 
class  of  correlations  as  the  timescale  correlations.  Notable  timescale  correlations  include  those  of  Ozawa17  and 
Plee18.  Ozawa’s  stability  correlation  is  actually  only  a  small  part  of  a  large  paper  on  the  properties  of  high-speed 
combustion  systems,  and  is  nicely  summarized  in  a  recent  paper  by  Driscoll19. 
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The  success  of  the  Damkohler  number  approach  depends  on  how  the  timescales  are  defined  and  calculated.  In 
Ozawa’s  correlation,  zytoli,  is  formed  from  the  inlet  velocity  and  the  characteristic  dimension  of  the  flameholder.  This 
timescale  is  said  to  represent  the  gas  residence  time  in  the  shear  layer.  The  chemical  timescale  is  described  as  the 
ignition  delay  time  of  the  mixture.  zchem  and  r/fell  are  used  to  form  a  Damkohler  number: 

Da  =  — (3) 

chem 

When  LBO  ®  is  plotted  against  Da,  a  typical  stability  loop  results.  Unlike  the  empirical  correlations,  Ozawa’s 
correlation  is  derived  using  multiple  data  sets  in  the  literature  for  different  conditions,  so  it  is  expected  to  apply  more 
generally  than  the  empirical  correlations.  Despite  the  improved  physical  description,  Ozawa’s  correlation  remains 
highly  empirical.  Ozawa’s  z chem  is  determined  by  plotting  zjimv  and  LBO  ®  data  and  working  backwards  to  determine 
zchem,  which  is  a  rather  indirect  method  for  determining  rctem. 

Plee  developed  one  of  the  most  advanced  timescale  correlations  in  the  literature18.  In  addition  to  chemical  and 
flow  timescales,  Plee’s  includes  timescales  to  account  for  mixedness  and  fuel  droplet  penetration  into  the  reaction 
zone.  For  well-mixed  conditions  with  a  fully  vaporized  fuel,  Plee’s  correlation  simplifies  to  a  form  very  similar  to 
Ozawa’s — stability  is  determined  based  on  competition  between  a  flow  and  chemical  timescale.  Like  Ozawa,  Plee 
uses  a  global  flow  velocity  and  flameholder  dimension  to  define  z/iow.  Plee  extends  the  development  of  the  chemical 
timescale  concept  beyond  Ozawa  by  identifying  that  ignition  delay  time  data  typically  follow  an  Arrhenius  form  and 
calculating  Tchem  directly.  Plee’s  approach  remains  highly  empirical  however,  in  that  the  Arrhenius  constants  and  the 
functional  relationship  between  Zjiow  and  zchem  used  to  predict  stability  are  empirically  derived. 

IV,  Limitations  of  Existing  Correlations 

Existing  correlations  are  unsatisfactory  for  several  reasons.  The  physical  justification  for  existing  correlations  is 
based  on  local  processes,  but  the  local  processes  are  estimated  from  global  parameters.  Calculating  local  parameters 
directly  would  increase  accuracy  and  make  the  methodology  more  general.  In  addition,  the  fundamental  assumption 
that  a  single  global  flow  and  chemical  timescale  determine  stability  is  an  oversimplification  for  modern 
flameholders.  For  a  multi-component  flameholder  with  a  complex,  three-dimensional  geometry,  multiple  timescales 
at  various  spatial  locations  may  be  important.  A  discussion  of  the  limitations  of  the  use  of  global  quantities  to 
predict  stability  has  been  discussed  by  Nair  and  Lieuwen4. 

Another  shortcoming  of  existing  stability  correlations  is  that  they  are  highly  dependent  on  empirical  data.  The 
empirically  derived  parameters  do  not  generally  extrapolate  well  to  new  configurations,  so  that  experimentation  is 
required  when  applying  an  existing  correlation  to  validate  or  modify  an  existing  correlation.  For  example.  King’s 
parameter  is  based  on  JP-4  fuel.  In  order  to  apply  this  correlation  to  a  configuration  involving  a  different  fuel,  such 
as  JP-10,  new  blowout  data  would  need  to  be  collected  in  order  to  empirically  derive  the  functional  relationship 
between  King’s  parameter,  Eq.  (2),  and  LBO  ®. 

A  final  limitation  is  that  the  conditions  for  the  experiments  that  form  the  basis  for  the  early  correlations  is  not 
always  well  defined.  There  is  often  confusion  between  inlet  velocity  and  the  mean  velocity  at  the  flameholder 
which  takes  into  account  the  blockage.  The  inlet  turbulence  field  is  often  not  provided  or  accurately  described,  as  is 
the  quality  of  the  mixedness  of  the  fuel  and  air.  The  lack  of  a  concise  description  of  the  conditions  for  which  the 
existing  correlations  were  derived  limits  the  ability  to  extrapolate  these  correlations  to  new  conditions.  Recent 
advances  in  diagnostics  and  CFD  have  greatly  improved  our  ability  to  resolve  with  high  spatial  and  temporal 
resolution  the  flow  fields  and  scalar  concentrations  in  these  types  of  flows. 

V,  Improved  Correlation  Methodology 

We  present  a  new  methodology  for  developing  stability  correlations.  The  same  fundamental  physical  description 
that  is  used  in  the  timescale  correlations  is  retained — stability  is  determined  using  a  Damkohler  number  formed  from 
a  relevant  flow  and  chemical  timescale.  Flowever,  unlike  existing  timescale  correlations,  the  timescales  are  not 
evaluated  using  global  quantities,  but  are  evaluated  locally.  In  addition,  instead  of  deriving  the  timescales 
empirically,  we  exploit  the  power  of  CFD  and  modern  chemical  kinetics  solvers  to  provide  timescale  data.  These 
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techniques  allow  for  evaluation  of  the  local  spatial  distribution  of  Da  in  the  reaction  zone,  as  opposed  to  using  a 
single,  global  Da  to  determine  stability.  The  computational  techniques  required  are  straightforward  and  have 
relatively  low  computational  cost. 

For  augmentors,  it  has  been  argued  that  the  relevant  flow  time  is  either  the  flow  residence  time  in  the  wake  or  the 
shear  layer  generated  by  the  flameholder.  In  our  approach,  the  turbulent  flow  field  in  the  wake  of  the  flameholder  is 
calculated  directly  using  unsteady  Reynolds-Averaged  Navier  Stokes  (RANS)  CFD  calculations  of  the  non-reacting 
flow.  For  a  particular  geometry,  this  procedure  can  be  repeated  over  a  range  of  Reynolds  numbers  to  generate  a 
correlation  for  flow  timescale  based  on  CFD  as  opposed  to  experimental  data. 

The  relevant  chemical  time  for  augmentor  stability  is  typically  described  as  an  ignition  delay  time.  Existing 
correlations  do  not  directly  calculate  or  estimate  this  timescale,  it  is  back-calculated  based  on  empirical  data.  In  our 
approach,  we  calculate  ignition  delay  times  computationally  using  reactor  models.  For  well-mixed  conditions,  the 
Perfectly-Stirred  Reactor  (PSR)  model  can  be  used.  For  poorly  mixed  conditions,  the  Partially-Stirred  Reactor 
(PaSR)  model  can  be  used.  Information  about  the  mixing  status  can  be  obtained  from  the  CFD  calculations.  Due  to 
the  low  dimensionality  of  reactor  models,  ignition  delay  data  can  be  calculated  relatively  easily,  even  for  complex 
chemical  mechanisms.  This  approach  allows  for  large  amounts  of  zchem  data  to  be  calculated  for  various  fuels, 
vitiation  levels,  and  inlet  pressures. 

To  summarize  our  approach,  low-order  CFD  calculations  are  used  to  determine  the  local  velocity  field  in  the 
wake  of  the  flameholder.  Reactor  models  are  used  to  calculate  chemical  timescales  for  the  flow.  A  Da  field  in  the 
wake  of  the  flameholder  is  calculated  and  used  to  determine  stability.  A  key  question  that  arises  when  evaluating 
stability  based  on  local  Da  values  is  what  location  or  location(s)  determine  stability?  This  is  an  open  question 
requiring  future  investigation.  Below  we  demonstrate  our  methodology  by  making  predictions  for  a  v-gutter 
flameholder  in  a  bluff-body  flameholder  rig. 


VI.  AFRL  In-House  Research  Rig 

We  demonstrate  the  methodology  by  generating  a  correlation  for  a  premixed,  propane-air,  v-gutter  configuration 
for  which  experimental  data  are  readily  available.  The  rig,  shown  in  Fig.  1  below,  has  a  rectangular  cross  section, 
5"  x  6",  so  that  the  flow  is  approximately  two-dimensional.  The  flameholder  considered  here  is  an  open  v-gutter 
with  1.5"  characteristic  dimension.  Our  collaborators  at  Wright-Patterson  Air  Force  Base  have  recently  performed 
lean  blowout  measurements  as  a  function  of  inlet  velocity  and  temperature  for  propane  fuel  for  this  rig.  Mean  inlet 
velocities  in  the  range  3.5  to  33  and  inlet  temperatures  in  the  range  300  to  530  K.  In  addition,  they  provided  inlet 
turbulence  intensity  data  and  have  recorded  a  large  volume  of  high-speed  image  data  for  observation  of  the  blowout 
dynamics.  An  example  is  superimposed  in  Fig.  1.  It  should  be  noted  that  although  we  assume  that  the  flow  is 
well-mixed  in  the  analysis  below,  there  is  some  uncertainty  as  to  whether  the  flow  was  actually  well-mixed  for  the 
data  presented  below.  These  results  should  therefore  be  treated  as  a  work  in  progress. 


Fuel  Injection 


Figure  1.  AFRL  in-house  research  rig  located  at  Wright-Patterson  Air  Force  Base. 


5 

American  Institute  of  Aeronautics  and  Astronautics 


VII.  Flow  Timescale 


We  calculated  the  flow  field  for  the  Augmentor  HIT  rig  using  unsteady,  two-dimensional  RANS  in  FLUENT® 
with  a  k-8  turbulence  model.  The  computational  domain  is  shown  in  Fig.  2  below.  Inlet  turbulence  intensity  was 
provided  experimentally,  although  the  turbulence  field  can  be  calculated  computationally  by  extending  the  domain 
further  upstream  in  the  flow.  FLUENT  simulations  were  performed  for  Re  in  the  range  10,000  to  55,000.  The 
structure  of  the  flameholder  wake  and  shear  layer  is  shown  in  Fig.  2  below,  which  shows  instantaneous  and  mean 
turbulence  intensity  (u  ’)  fields  for  Re  =  40,000.  CFD  calculations  were  validated  based  on  the  Strouhal  number  of 
the  characteristic  vortex  shedding  (St  =  0.275)  and  by  comparing  mean  and  RMS  velocity  profiles  with  LDV  data. 
In  general,  agreement  between  our  CFD  calculations  and  experimental  data  for  the  flow  field  is  good. 


Figure  2.  Instantaneous  (top)  and  mean  u  ’  fields  for  the  Augmentor  HIT  rig,  Re  =  40,000,  showing  the  vortex 
shedding  structure  typical  of  a  flameholder  wake. 

For  the  current  demonstration,  we  used  the  ratio  of  turbulent  kinetic  energy  (k)  and  dissipation  (e)  to  form  the 
flow  timescale,  Zjjow  =  k/e.  This  timescale  is  representative  of  the  local  rate  of  turbulent  mixing.  It  should  be  noted 
that  this  definition  of  r//o„,  has  not  been  rigorously  validated  but  is  used  here  as  a  reasonable  starting  point  for 
demonstrating  the  overall  methodology.  The  k/e  field  for  the  shear  layer  formed  at  the  downstream  lip  of  the 
flameholder  is  shown  in  Fig.  3  below,  k/e  values  reach  a  minimum  in  the  shear  layer  close  to  the  lip  of  the 
flameholder.  For  this  simple  rig  configuration,  we  assume  that  this  region  controls  stability  and  we  evaluate  Zfiowfo r 
Da  in  this  critical  region.  For  the  conditions  considered  here  zjiow  is  in  the  range  0.2  to  2  ms. 

In  the  current  demonstration,  non-reacting  CFD  calculations  are  used  to  calculate  zjjow.  The  error  associated  with 
this  assumption  is  currently  unknown.  In  the  future  we  plan  to  quantify  this  error  by  comparing  zjiow  values 
generated  using  non-reacting  and  simple  reacting-flow  CFD  calculations. 


Critical 

Region 


Figure  3.  k/e  field  close  to  the  trailing  edge  of  the  flameholder.  k/s  values  are  minimum  in  the  shear  layer 
close  to  the  flameholder  and  reaching  values  ~0.5  ms. 
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VIII.  Chemical  Timescale 

t chem  values  are  calculated  based  on  ignition  delay  time  calculations  made  using  using  a  perfectly  stirred  reactor 
model.  Calculations  are  made  in  CHEMKIN 8  using  the  San  Diego  propane  mechanism.  All  calculations  are  for 
propane-air  at  1  atm.  Ignition  delay  calculations  are  highly  sensitive  to  initial  temperature.  Plee18  has  argued  that 
the  maximum  temperature  in  the  system  should  be  used  as  the  initial  temperature,  although  how  to  determine  this 
temperature  is  another  open  question.  Here  we  use  adiabatic  flame  temperature  ( Tad)  as  the  initial  temperature  in 
order  to  be  as  objective  as  possible.  Using  Tad  results  in  the  fastest  possible  values  for  zchem  and  probably  under¬ 
predicts  ichem  relative  to  actual  conditions. 

An  ignition  delay  time  curve  is  shown  for  propane  air  for  ®  =  0.3  in  Fig.  4  below.  Ignition  delay  is  calculated 
based  on  peak  OH  concentration,  which  is  a  common  convention.  For  this  very  lean  equivalence  ratio,  zchem  is 
relatively  slow,  ~2  ms.  zch(!m  is  calculated  as  a  function  of  ®  in  this  way  to  generate  a  curve  relating  zchem  and  ®  for 
lean  propane-air.  Fig.  5  below.  The  relationship  is  highly  nonlinear. 

Propane/Air  Autoignition,  cp  =  0.3 
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Figure  4.  Propane-air  ignition  delay  calculation,  ®  =  0.3,  T„  =  Tad. 
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Figure  5.  rc/,c,„  as  a  function  of  ®  for  lean  propane-air. 
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IX.  LBO  Predictions 


For  the  AFRL  in-house  rig  configuration,  we  evaluate  Da  in  the  critical  region  as  shown  in  Fig.  3.  Our  stability 
limit  is  simply  Da  =  1 .  Our  predictions  are  plotted  against  measured  values  and  values  generated  using  the  existing 
correlations  described  above  in  Fig.  6  below.  Values  for  King  are  not  shown  because  King’s  correlation  does  not 
apply  to  the  temperature  range  for  the  FIIT  rig  data  and  returns  non-physical  (rich)  LBO  <t>.  As  expected,  the 
timescale  correlations  of  Ozawa  and  Plee  outperform  the  empirical  correlations  of  King  and  DeZubay.  The  existing 
timescale  correlations  are  significantly  offset  in  ®  from  the  measured  values.  The  stratification  in  the  Plee  data  is 
grouped  by  inlet  temperature. 

The  current  method  does  a  good  job  of  predicting  this  configuration.  This  result  should  be  treated  with  caution 
due  to  the  potential  mixing  issue  described  above.  The  reported  LBO  ®  values  are  unusually  low,  with  many  data 
points  falling  below  ®  =  0.4,  lending  credence  to  the  possibility  of  a  fuel  mixing  problem.  In  addition,  many 
assumptions  were  made  in  this  demonstration  of  our  methodology,  such  as  perfect  mixing  and  assuming  an  initial 
temperature  of  Taj  for  ignition  delay  calculations,  that  could  prove  faulty  when  this  methodology  is  applied  to  more 
realistic  configurations.  We  would  like  to  emphasize  that  we  are  presenting  an  improved  methodology,  and  that 
there  are  many  open  questions  regarding  how  the  methodology  should  be  applied  for  a  particular  configuration.  For 
the  current  configuration,  our  methodology  was  developed  and  implemented  so  that  LBO  predictions  would  match 
known  experimental  values.  A  useful  design  tool  must  be  able  to  accurately  predict  various  configurations  without 
a  priori  knowledge  of  the  outcome.  It  is  also  very  interesting  to  note  that  the  existing  correlations  predict  LBO  over 
a  broad  range  of  values,  indicating  the  predictive  inconsistency  of  the  existing  approaches  independent  of  our 
predictions. 
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Figure  6.  Lean  blowout  ®  predictions  using  the  current  method  and  several  existing  correlations  for  the 
AFRL  in-house  rig. 


In  this  demonstration,  stability  was  determined  based  on  Da  values  at  a  particular  location  in  the  flow.  For 
future  analyses  of  more  complex  geometries,  we  plan  to  develop  a  more  advanced  approach  where  the  local  spatial 
distribution  of  Da  is  used  to  determine  stability,  as  opposed  to  a  single  local  value  for  Da.  In  this  way  more 
advanced  criteria  can  be  applied  to  determine  stability.  When  implementing  this  approach,  we  plan  to  form  a 
probability  density  function  (PDF)  for  Da  describing  the  spatial  distribution  of  Da  in  the  reaction  zone.  Stability 
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predictions  could  then  be  made  based  on  this  PDF.  For  example,  if  local  extinction  were  to  occur  for  a  certain 
volume  of  the  volume  of  the  flow  (i.e.,  if  a  certain  volume  has  Da  <  1),  blowout  would  occur. 

Similarly,  a  model  incorporating  unsteady  flow  effects  could  be  developed  based  on  both  the  spatial  and 
temporal  variations  in  Da  in  the  reaction  zone.  For  flames  where  vorticity  plays  a  key  role  in  blowout,  temporal 
variations  in  Da  may  be  as  critical  as  spatial  variations  in  determining  stability. 

X.  Challenges  and  Future  Work 

Progress  in  this  research  program  to  date  reflects  the  scope  of  a  relatively  small  Phase  I  SBIR  research  program. 
The  methodology  has  been  developed  and  demonstrated  for  a  single,  simplified  configuration.  This  exercise  has 
raised  many  important  questions  on  how  this  methodology  should  be  implemented.  What  are  the  critical  location(s) 
in  an  augmentor  that  determine  stability?  What  temperature  should  be  used  when  making  ignition  delay 
calculations?  Can  cold-flow  CFD  simulations  be  used  to  calculate  Tjlow  values  with  reasonable  accuracy?  When  can 
the  perfect  mixing  assumption  be  applied,  and  when  should  a  PaSR  reactor  be  used?  Flow  should  fuel  penetration 
effects  be  incorporated  into  the  model?  In  the  future  we  hope  to  address  these  questions  by  implementing  our 
approach  on  a  broad  range  of  configurations,  including  other  rig  configurations  and  realistic  combustor 
configurations. 
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